Electrical Conductivity of Lithium at Mbar Pressures 
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We report measurements of the electrical conductivity of a liquid alkali metal - lithium - at pressures 
up to 1.8 Mbar and 4-fold compression, achieved through shock compression experiments. We find 
that the results are consistent with a departure of the electronic properties of lithium from the 
nearly free electron approximation at high pressures. 



PACS numbers: 72.15.-v, 62.50.+p, 71.22. +i 

At one atmosphere lithium is the archetype of a "sim- 
ple" metal, i.e. one in which the electronic valence states 
are well separated energetically from the tightly bound 
core states. It has been known for a long time that 
its electronic properties at these ambient conditions are 
well described within the nearly free electron picture |L| , 
which is a cornerstone of the theory of simple metals. 
While it has been generally assumed that increased pres- 
sure would only improve the accuracy of this descrip- 
tion, in recent years theoretical calculations suggested 
that the opposite may be true M. This has culminated 
very recently with the striking prediction that lithium, 
long viewed as the simplest of all metals, exhibits at Mbar 
pressures a paired ground state with a semiconducting or 
insulating character 0]. The ensuing experimental work 
has provided results on the decrease of optical reflectiv- 
ity upon compression 0,p| , on the increase of the electri- 
cal resistivity with pressure up to 0.6 Mbar, and on 
the existence of symmetry-breaking structural transitions 
around 0.5 Mbar m, all in broad agreement with this 
prediction. It is perhaps worth noting that even before 
the nearly free electron model was proposed, Bridgman's 
high pressure experimental work showed that the resistiv- 
ity of lithium increases with pressure, both in the solid 
and in the liquid [pi. Here we report measurements oi 
the electrical resistivity of lithium at pressures up to 1.8 
Mbar, achieved through shock compression experiments. 
The results are consistent with a departure of the elec- 
tronic properties of lithium from the nearly free electron 
approximation at high pressures, and with ionic pairing 
correlations in the Mbar regime (IMbar — lOQGPa). 

We measured the electrical resistance of high purity 
99.995% lithium samples quasi-isentropically compressed 
starting from the solid with density do = .531g / cm 3 at 
room temperature and Pq = lbar. The quasi-isentropic 
compression was achieved through multiple reflections 
of a shock wave M between two sapphire single-crystals 
which encapsulate the lithium sample. The initial shock 
wave was generated by the impact of Al or Cu projectiles 
moving at 3 to 7km/ s onto the experimental cell (target) 
containing the samples. The targets were specifically de- 
signed to create and maintain steady state conditions at 
the final pressure for time durations exceeding 100ns, 
during which the measurements were taken. 

Due to the extreme reactivity of lithium, the exper- 



imental cells were assembled inside an Ar atmosphere 
glove-box, and designed to be hermetically sealed after 
assembly. Pre-shock sample resistances were within 5% 
of the ideal room temperature values and were monitored 
continuously after the targets were removed from the con- 
trolled environment to ensure that the sample quality did 
not degrade. In order to increase the electrical resistance 
and better control current flow, we used parallelipipedic 
samples with dimensions length x width = 20mm x 3mm 
and thicknesses varying between .2mm and .3mm. The 
10 : 1 and bigger aspect ratios used insured that lateral 
rarefaction effects were minimal. 




FIG. 1. Schematic representation of the experimental cell 
(not to scale). (1) Al target body; (2) single-crystal sapphire 
shock anvils; (3) Lithium sample; (4) Polyethilene filler - pro- 
vides the electrical insulation around the sample; (5) Al or 
Cu foil matching the projectile material to provide screening 
of electrical discharges produced upon impact between the 
metallic flyer and the target; (6) triggering pins; (7) 4 Cu 
electrodes for the conductivity measurements. 

The electrical resistance of the lithium samples was 
measured using a four probe technique, which virtually 
eliminates the need for the contact resistance correction. 
The electrodes were gold-plated, oxygen-free Cu wires, 
inserted through the back sapphire anvil, see Fig.[j]. The 
diameter and positioning of the outer electrodes were op- 
timized in order to provide uniform current injection into 
the lithium sample, and to minimize the contact resis- 



tance and distortions of the pressure profile in the mea- 
surement region. The voltage probes - inner electrodes - 
were very thin and placed 3 — 6mm apart as needed to 
maximize the sample resistance and the accuracy of the 
voltage measurements. Special attention was given to the 
soldering process and shielding of the electrodes in order 
to improve the signal-to-noise ratios. Triggering of the 
data acquisition system was provided by time-of-arrival 
shock sensors placed outside the sample space in order to 
eliminate any interference with the measurements. The 
details of the electronic circuitry used in these experi- 
ments are similar with the ones described in ||l0||. 
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FIG. 2. Example of a typical experimental trace. Measured 
sample resistance (scale on left). vs. time - continuous line; cal- 
culated pressure (scale on right). vs. time - dashed line. 

The sample pressures at final, steady-state conditions 
were determined with 1% accuracy from the measured 
projectile velocity using the shock impedance matching 
technique flilj . The time dependence of the pressure dur- 
ing compression was calculated using a one-dimensional 
hydrodynamic code in which the projectile and shock- 
anvils were modeled by Mie-Gruneisen equations of state 
(EOS) @. The EOS for the lithium sample was ob- 
tained by least-squares fits of a ratio of polynomials to 
tabular data based on extensive Shockwave and static 
compression experimental results |l3]. We note that all 
the conductivity results reported in this paper are in the 
liquid region of the phase diagram. The densities and 
temperatures attained in the experiments range approx- 
imately from 2 to 4 times the normal density, and from 
2000 A' to 7000 A. 

In the experiment wc monitored the voltage drop 
across the lithium sample as the pressure increased and a 
known constant current was passed through, see Fig. |2j. 
The accuracy in the current and voltage measurements 
was higher than 1% and 5%, respectively. The lithium 
electrical resistivity was determined from the sample re- 
sistance using R = px l/S, where the distance I between 
the voltage probes was measured to 0.1% precision. The 



cross-section of the sample, S — w x th, was determined 
by the known width w and calculated thickness th, with a 
conservatively estimated precision of 10% arising mainly 
from density uncertainties. 
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FIG. 3. Measured electrical conductivity of lithium as a 
function of pressure. Open symbols - data extracted from the 
final, steady-state conditions; crossed circles - data extracted 
from intermediate states. Note: IMbar — WOGPa. 

The measured electrical resistivity exhibits three main 
regimes as a function of pressure, see Fig. [| Up to 
approximately IMbar and 3 — fold compression the elec- 
trical conductivity of lithium decreases steeply from its 
10 5 (£lcm)~ 1 value at normal pressure to about 3.6 x 
10 3 (ficm) -1 at IMbar, which is consistent with increased 
scattering of the charge carriers in the compressed ma- 
terial. This regime is followed by a much slower vari- 
ation of the conductivity with pressure, which extends 
to ~ l.QMbar and almost 4 — fold compression. At 
even higher pressures, although the error bars overlap, 
the electrical conductivity appears to depart from the 
previous slow varying behavior. 

In order to interpret the results we turn to the general 
framework provided by the electron-ion pseudopotential 
formalism and the Ziman conductivity theory p4| . The 
pseudopotential concept takes into account the near can- 
cellation of the Coulombic interaction between the va- 
lence electrons and ions inside the core, due to orthog- 
onality and exclusion effects. While generally non-local 
and energy dependent, in its most intuitive representa- 
tion it is just a local potential |l5| . Under the assumption 
of a weak pseudopotential the linear response theory ap- 
plied to the uniform electron gas yields the total energy of 
the system and defines effective ion-ion and electron-ion 
interactions p6| , p7| . 

The effective ion-ion pair potential is the sum of a di- 
rect Coulombic interaction, and an indirect contribution 
due to the polarization of the electrons, 
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and the screened electron-ion potential is, 
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where Z is the number of valence electrons per atom, x{k) 
is the static response function, w(k) is the bare electron- 
ion pseudopotential, Xo{k) is the Lindhard polarizability, 
e(k) is the static dielectric function and G{k) is the local 
field factor of the electron gas, accounting for the ex- 
change and correlation effects between the electrons [[L7| . 
The electrical resistivity follows from the assumption 
of a degenerate electron gas and weak electron-ion scat- 
tering, yielding the Ziman formula |l7|-|lS|l : 
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The deviations of the electronic density of states from 
the free electron values, due to the pseudopotential, are 
accounted for through the effective mass m* po|~f22[; S 
is the liquid-structure factor, y = k/2kp, kp is the Fermi 
wavevector and ao is the Bohr radius. 

As remarked in p] the strongly non-local character of 
the pseudopotential generates significant deviations of 
the electronic structure from the free electron form at 
high densities, reflected in part in the decrease of the oc- 
cupied bandwidth, i.e. increase of the effective mass. 
We note that, under suitable assumptions, Eq. @ al- 
lows us to estimate the effective mass m* , by comparing 
with the measured experimental results. For the sake 
of simplicity and to gain additional insight (see below), 
we use the local empty core potential of Ashcroft p5| , 
w(r) — for r < R c and —Ze 2 /r above R c . However, 
we do not expect that the use of a non-local potential 
for the calculation of the scattering integral will quali- 
tatively alter the results |2lf|. Values of R c typical for 
lithium range from 1.06a. u. |18| to 1.44a.u. p3fl , depend- 
ing on the local field correction employed. We use the 
homogeneous electron gas local field factor G(k) deter- 
mined by diffusion Monte Carlo simulations [E4[ , and find 
that R c = 1.26a. it. reproduces the experimental values of 
the conductivity of liquid lithium close to its triple point 
[ 25 1 with an effective mass only slightly bigger than unity 
[21|. We also note that the normal pressure effective ion- 
ion potential so obtained compares favorably with other 
potentials that have been used to model liquid lithium 
p3| . The liquid-structure factor S(k) of liquid metals 
has been usually assumed to be well reproduced by the 
hard-sphere model with some appropriate packing frac- 
tion |lq]. Here we determine S(k) numerically based on 



the effective ion-ion interaction $(r), using the pertur- 
bative hypernetted-chain equation (PHNC) {||. PHNC 
has been successfully applied to various model systems, 
including liquid metals p7[ , and should be particularly 
appropriate at high densities. We note here that the 
temperature dependence of the scattering integral is em- 
bedded in S(k). We calculate for example that a 20% 
increase in temperature yields a 4% decrease in the scat- 
tering integral at 2-fold compression, and only about 2% 
at 4-fold compression. This is due to the fact that at 
high densities the liquid structure is determined mainly 
by the ion- ion repulsions |2q ], whose role only increases 
with increased compression. 

The results of the effective mass calculations based on 
the Ziman theory with R c = 1.26a.u. are summarized in 
Fig. U We note that the initial steep decrease in con- 
ductivity with pressure corresponds to an increase of m* 
to values similar with the ones found in [f|, as lithium 
becomes less free-electron like upon compression. In this 
pressure range the reduction in conductivity is driven 
mostly by increased scattering due to decreased inter- 
atomic spacing, although the core electronic states ap- 
pear to already come into play as evidenced by changes 




FIG. 4. Electronic effective mass as a function of r s 
(r s — (3/47ro e ao) , where n e is the valence electrons num- 
ber density, with one electron/atom; m is the electron mass). 
Open circles - estimated effective mass for R c = 1.26a. it.; 
dashed line - extrapolation of the lower pressures regime; tri- 
angle - effective mass reported in Ref. [3]. Inset shows the 
variation of the effective core radius with r s along the dashed 
line. The units for R c and r s are a.u.. 

in the effective mass. In the next regime the slow vari- 
ation of the conductivity, see Fig. g, is accompanied by 
a decrease of the Ziman estimated effective mass, Fig. 
H. It is very likely in fact that the actual effective mass 
]22| increases monotonically under compression past the 
predicted pairing instability [pf . This is the behavior ob- 
served in LDA calculations of a dense lithium monolayer, 
which is believed to behave similarly with bulk lithium 
|E9| . We note therefore that the slow variation of the con- 



ductivity, along with the decrease of the Ziman estimated 
effective mass, point to a shrinking ionic core, which bal- 
ances the increase in scattering due to increased density. 
In the final regime, as the conductivity appears to drop, 
the effective mass increases again. In the framework of 
an empty core model this suggests that the exclusion- 
ary effects of the ions on the scattering of the valence 
electrons become dominant over the density driven core 
decrease. To make these observations more quantitative, 
we extrapolate the lower density effective mass to the 
high density value reported in pfl , and determine the core 
radius R c that reproduces this trend - see inset to Fig. 

1 

As discussed in y,Q in connection to several alkali 

metals, the importance of the core electronic states in- 
creases at high pressures, and along with it the magni- 
tude of the pseudopotential. As a result, in the solid, 
symmetry-breaking distortions leading to ionic pairing 
may become energetically favorable |3|]. In the liquid, a 
rising pseudopotential, i.e. a decreasing effective ionic 
core, see Fig. ||, should also mediate a strengthening of 
pairing correlations. Due to the strong exclusionary ef- 
fects of the ions on the valence electrons at high densities, 
such correlations maximize in turn the spaces available 
for the electrons in the interstitial regions, away from 
the areas of strong core overlap in-between the ions, with 
beneficial effects on the energy. This exclusion of the va- 
lence electrons by the ions seems to be consistent with 
the behavior of the conductivity at high compressions. 
Given the expected small effect of the temperature on the 
conductivity at high densities, the apparent conductivity 
drop and the behavior of the ionic core at the highest 
pressures could be interpreted as a decrease of the over- 
all volume available for the electrons. 

It may be interesting to see if the exclusionary effects 
mentioned above, that lead to a very non-uniform distri- 
bution of the valence electrons, translate into the analog 
of classical depletion forces [[Hj , enhancing the ionic pair- 
ing correlations in a mixture, e.g. LiH. Higher tempera- 
tures would ultimately destroy pairing correlations in the 
liquid. Estimates of the temperatures required should be 
possible, but need to rely on more detailed calculations, 
that could be tested against the experimental results pre- 
sented here. 

We thank N. Ashcroft, J. Neaton, G. Galli and R. 
Cauble for useful discussions, and D. Young for the 
lithium EOS table. M.B. warmly thanks A. Mitchell 
for sharing his life-long experience on Shockwave exper- 
iments. We gratefully acknowledge the LDRD office for 
financial support. We also thank S. Caldwell, W.P. Hall, 
N.Hinsey, K. Stickle, L. Raper and T. Uphaus for assis- 
tance at the gas- gun facility, and L. Walkley for help with 
the set-up of the glove-box. This work was performed 
under the auspices of the U. S. Department of Energy 
by University of California Lawrence Livermore National 
Laboratory under Contract No. W-7405-Eng-48. 



' To whom correspondence should be addresed: 

basteal@llnl.gov 
[1] F. Seitz, Phys. Rev. 47, 400 (1935). 
[2] J.C. Boettger, S.B. Trickey, Phys. Rev. B 32, 3391 

(1985). 
[3] J.B. Neaton, N.W. Ashcroft, Nature 400, 141 (1999). 
[4] Y. Mori, C.-S. Zha, A. Ruoff, in Science and Technology 

of High Pressure, Proceedings of AIRAPT-17, edited by 

M.H. Manghnani, W.J. Nellis and M.F. Nicol, p. 421, 

Universities Press, Hyderabad, India (2000). 
[5] V.V. Struzhkin, R.J. Hemley, H.K. Mao, Bull. Am. Phys. 

Soc. 44, 1489 (1999). 
[6] V.E. Fortov, V.V. Yakushev, K.L. Kagan, I.V. 

Lomonosov, V.I. Postnov, T.I. Yakusheva, JETP Lett. 

70, 628 (1999). 
[7] M. Hanfland, K. Syassen, N.E. Christensen, D.L. 

Novikov, Nature 408, 174 (2000). 
[8] P.W. Bridgman, Phys. Rev. 17, 161 (1921). 
[9] M. Bastea, A.C. Mitchell, W.J. Nellis, Phys. Rev. Lett. 

86, 3108 (2001). 
[10] A.C. Mitchell, R.N. Keeler, Rev. Sci. Instr. 39, 513 

(1968). 
[11] A.C. Mitchell, W.J. Nellis, Rev. Sci. Instr. 52, 347 

(1981). 
[12] S.P.Marsh, LASL Shock Hugoniot Data, University of 

California Press, Berkeley (1980). 
[13] DA. Young, Phase Diagram of the Elements, Univer- 
sity of California Press, Berkeley (1991), and references 

therein. 
[14] J.M. Ziman, Phil. Mag. 6, 1013 (1961). 
[15] N.W. Ashcroft, Phys. Lett. 23, 48 (1966). 
[16] N.W. Ashcroft, D.C. Langreth, Phys. Rev. 159, 500 

(1967). 
[17] M. Shimoji, Liguid Metals, Academic Press, London 

(1977). 
[18] N.W. Ashcroft, J. Phys. C 1, 232 (1968). 
[19] D.J. Stevenson, N.W. Ashcroft, Phys. Rev. A 9, 782 

(1974). 
[20] M. Itoh, M. Watabe, J. Phys. F: Met. Phys. 14, L9 

(1984). 
[21] A.B. Van Oosten, W. Geertsma, Physica 133B, 55 

(1985). 
[22] T. Chan, L.E. Ballentine, Can. J. Phys. 50, 813 (1972). 
[23] M. Canales, J. A. Padro, Phys. Rev. E 56, 1759 (1997). 
[24] S. Moroni, D.M. Ceperley, G. Senatore, Phys. Rev. Lett. 

75, 689 (1995). 
[25] T.C. Chi, J. Phys. Chem. Ref. Data 8, 339 (1979). 
[26] H.S. Rang, F.H. Ree, J. Chem. Phys. 103, 3629 (1995). 
[27] H.S. Rang, Phys. Rev. B 60, 6362 (1999). 
[28] N.W. Ashcroft, J. Lekner, Phys. Rev. 145, 83 (1966). 
[29] A. Bergara, J.B. Neaton, N.W. Ashcroft, Phys. Rev. B 

62, 8494 (2000). 
[30] J.B. Neaton, N.W. Ashcroft, Phys. Rev. Lett. 86, 2830 

(2001). 
[31] see, for example, V. Vlachy, Langmuir 12, 2881 (1996). 



